Stabilization of Hydrodynamic Flows by Small Viscosity Variations 
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Motivated by the large effect of turbulent drag reduction by minute concentrations of polymers 
we study the effects of a weakly space-dependent viscosity on the stability of hydrodynamic flows. 
In a recent Letter [Phys. Rev. Lett. 87, 174501, (2001)] we exposed the crucial role played by a 
localized region where the energy of fluctuations is produced by interactions with the mean flow (the 
"critical layer" ) . We showed that a layer of weakly space-dependent viscosity placed near the critical 
layer can have a very large stabilizing effect on hydrodynamic fluctuations, retarding significantly 
the onset of turbulence. In this paper we extend these observation in two directions: first we show 
that the strong stabilization of the primary instability is also obtained when the viscosity profile 
is realistic (inferred from simulations of turbulent flows with a small concentration of polymers). 
Second, we analyze the secondary instability (around the time-dependent primary instability) and 
find similar strong stabilization. Since the secondary instability develops around a time-dependent 
solution and is three-dimensional, this brings us closer to the turbulent case. We reiterate that the 
large effect is not due to a modified dissipation (as is assumed in some theories of drag reduction), 
but due to reduced energy intake from the mean flow to the fluctuations. We propose that similar 
physics act in turbulent drag reduction. 
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I. INTRODUCTION 

This paper is motivated by the dramatic effects that 
are observed with the addition of small amounts of poly- 
mers to hydrodynamic flows. While interesting effects 
were discussed in the context of the transition to turbu- 
lence, vortex formation and turbulent transport Q, the 
phenomenon that attracted the most attention was, for 
obvious reasons, the reduction of friction drag by up to 
80% when very small concentrations of long-chain poly- 
mers were added to turbulent flows j| [j| . In spite of the 
fact that the phenomenon is robust and the effect huge, 
there exists no accepted theory that can claim quantita- 
tive agreement with the experimental facts. Moreover, 
it appears that there is no mechanistic explanation. In 
the current theory that is due to de Gennes [Q, ^| one ex- 
pects the Kolmogorov cascade to be terminated at scales 
larger than Kolmogorov scale, leading somehow to an in- 
creased buffer layer thickness and reduced drag, but how 
this happens and what is the fate of the turbulent energy 
is not being made clear. 

In a recent Letter Q we proposed that the crucial issue 
is in the production of energy of hydrodynamic fluctua- 
tions by their interaction with the mean flow. For the 
sake of concreteness we examined a Poiseuille laminar 
flow and its loss of linear stability, and showed how small 
viscosity contrasts lead to an order of magnitude retar- 
dation in the onset of instability of "dangerous" distur- 
bances. Specifically, we considered a flow in a channel of 
dimensionless width 2, in which there are two fluids: one 
fluid of viscosity [i\ flows near the walls, and the other 
fluid of viscosity /i2 flows at the center, see Fig. |l|. The 
viscosities differ slightly, for example we considered (in 
dimensionless units) /12 = 1 and m = ^1/^2 = 0.9. The 
main ingredient of the calculation was that all the vis- 
cosity difference of 0.1 concentrated in a "mixed" layer 
of width 0.10. The motivation behind these numbers 
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FIG. 1: Schematic of the flow: the fluid near the walls has a 
viscosity (ii, and that flowing at the center is of viscosity fj,2. 
In the mixed layer (of width g) the viscosity varies gradually 
between pi and fj.2- The parameter p controls the position of 
the mixed layer. For simplicity we neglect the down-stream 
growth in q. 



was the observation that the inferred effective viscosity 
in polymer drag reduction increases towards the center 
by about 30% over about a 1/3 of the half-channel [Q. 
With our choice we have comparable viscosity gradients 
in the mixed layer. 

In this model everything was explicitly calculable. The 
main point of our analysis (see Sect. |n] for further details) 
was that there exists a position in the channel where the 
velocity of the mean flow is the same as the velocity of 
the most dangerous primary instability. Below wc refer 
to the layer around this position as the "critical layer". 
If we placed the mixed layer in the vicinity of the critical 
layer, we got a giant effect of stabilization. Analyzing 
this phenomenon, we demonstrated that nothing special 
happened to the dissipation. Rather, it was the energy 
intake from the mean flow to the unstable mode that was 
dramatically reduced, giving rise to a large effect for a 
small cause. In this paper we extend these observations in 
two directions. In Sect. ||, after reviewing the results of 
the simple model, we extend the analysis of the primary 
instability to a case in which the viscosity profile is that 
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inferred from direct numerical simulations of turbulent 
channel flow of a dilute polymeric solution 0. We will 
see that very similar effects are found. In other words, 
one does not need to put by hand the region of viscosity 
variation in the vicinity of the critical layer. When we 
have a continuous variation of the viscosity in the region 
near the wall, the effect is the same, since it is only crucial 
that there will be some space dependence of the viscosity 
in the critical layer, which is usually not too far from the 
wall. 

A possible criticism of our results can be that the pri- 
mary instability is still too far from typical turbulent 
fluctuations. This is in particular true since the most 
unstable primary modes are 2-dimensional, whereas typ- 
ical turbulent fluctuations are 3-dimensional. For these 
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reasons we present in Sect. Ill the analysis of the effect 



of small viscosity variations on the secondary instability, 
for which the most "dangerous" modes are 3-dimensional. 
The tactics are similar to those taken for the primary in- 
stability. First we discuss the effects of a mixed layer put 
at the "right" place in the channel, and second we show 
that continuous viscosity profiles do exactly the same. 
We find again the giant effect of stabilization for rela- 
tively small viscosity variations, lending further support 
to our proposition that similar effects may very well play 
a crucial role in turbulent drag reduction. In Sect. |\J we 
present concluding remarks and suggestions for the road 
ahead. 



II. PRIMARY INSTABILITY OF POISEUILLE 
FLOW 

It is well known that parallel Poiseuille flow loses linear 
stability at some threshold Reynolds number Re=Reth 
(close to 5772). It is also well known that the instability 
is "convective" , with the most unstable mode having a 
phase velocity c p . Analytically it has the form 

4> p (x,y,t) = ^{(j) p (y) exp[ik p (x - c p t)} + c.c.} exp(^pt) , 

.(1) 

where the subscript p stands for the primary instability, 
<f>(x, y, t) is the disturbance streamfunction and 4>(y) is 
the complex envelope of <j>(x, y, t). We have chosen x and 
y as the streamwise and wall-normal coordinates respec- 
tively, k as the streamwise wavenumber of the distur- 
bance and t as time. j p is the growth rate of the primary 
instability. What is not usually emphasized is that the 
main interactions leading to the loss of stability occur in 
a sharply defined region in the channel, i.e. at the crit- 
ical layer whose distance from the wall is such that the 
phase velocity c is identical to the velocity of the mean 
flow somewhere within this layer. It is thus worthwhile 
to examine the effect on the stability of Poiseuille flow of 
a viscosity gradient placed in the vicinity of the critical 
layer. This will provide us with a very sharp understand- 
ing of the mechanism of the stabilization of the flow by 
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FIG. 2: The dependence of the threshold Reynolds number 
Reth on the position of the viscosity stratified layer for m = 
0.9. The dashed line pertains to the neat fluid. Note the huge 
increase in Kth within a small range. This occurs when the 
stratified layer overlaps the critical layer. 



viscosity variations. In the following subsection we will 
examine the case of continuous viscosity profiles. 



A. Mixed Layer 

A report of the results of this subsection was provided 
in ||. We examine a channel flow of two fluids with 
different viscosities /ii and /i2, see Fig. 0. 

The observation that we want to focus on is shown in 
Fig. 0: the threshold Reynolds number for the loss of 
stability of the mode as in Eq. (|l|) depends crucially on 
the position of the mixed layer. When the latter hits the 
critical layer, the threshold Reynolds number for the loss 
of stability reaches as much as 88000. In other words, 
one can increase the threshold of instability for a given 
mode 15 times, and by making the mixed layer thinner 
one can reach even higher threshold Reynolds values. In 
H we analyzed the physical origin of this huge sensitivity 
of the flow stability to the profile of the viscosity. 

The stability of this flow is governed by the modified 
Orr-Sommerfeld equation S 
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k^p) (U-c p - i 7p ) - U"cj> p ] 

^( 4 ) + 2^;" + (//'-2fc»< 



(2) 



in which U(y) is the basic laminar velocity, and fj, 
is a function of y. The boundary conditions are 
4> P {±1) = 4>' p (±l) — 0. All quantities have been non- 
dimensionalised using the half-width H of the channel 
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FIG. 3: Profiles of the normalized viscosity fi(y) and nor- 
malized velocity U(y) and the second derivative U (y) for 
m = 0.9 (solid lines) and m = 1.0 (dashed lines). The mixed 
layer is between the vertical dashed lines. 



It can be seen, comparing the mean profile U(y) to 
that of the neat fluid (cf. Fig. |J), that nothing dra- 
matic happens to this profile even when the mixed layer 
is chosen to overlap a typical critical layer. Accord- 
ingly we need to look for the origin of the large effect 
of Fig. 2 in the energetics of the disturbances. To do 
so, recall that the streamwise and normal components 
of the disturbance velocity u p (x, y, t) and v p (x, y, i) may 
be expressed via streamfunction as usual: U p (x, y, t) = 
d(f)p/dy, and v p (x,y,t) = —d(f> p /dx. These functions 
may be written in terms of complex envelopes similar to 
Eq. (|) : 

u p (x,y,t) = i{u p (?/)exp[ifcp(a;-Cpi)]+c.c.}exp(7pt)(9) 
v p (x,y,t) = i{u p (y)exp [ik p (x - c p t)] + c.c.} exp(M) . 

The pressure disturbance p p is defined similarly. 

Define now a disturbance of the density of the kinetic 
ene rgy of the primary instability 



and the centerline velocity Uq as the length and velocity 
scales respectively. The Reynolds number is defined as 
Re= pUoH/ /j,2, where p is the density (equal for the two 
fluids). The primes stand for derivative with respect to y. 
At y = 0, we use the even symmetry conditions <fi(0) = 1 
and (f>'(0) = 0, as the even mode is always more unstable 
than the odd. 

Since the flow is symmetric with respect to the channel 
centreline, we restrict our attention to the upper half- 
channel. Fluid 2 occupies the region < y < p. Fluid 1 
lies between p + q < y < 1 . The region p < y < p + q 
contains mixed fluid. The viscosity is described by a 
steady function of y, scaled by the inner fluid viscosity 

P2- 

p(y) = 1, for 0<y< P , (3) 

fi(y) = l + (m-lK 3 [lO-15e + 6£ 2 ] , 0<£<1,(4) 

fi(y) = m , for p + q < y < 1. (5) 

Here £ = {y ~ p)/q is the mixed layer coordinate. We 
have assumed a 5th-order polynomial profile for the vis- 
cosity in the mixed layer, whose coefficients maintain the 
viscosity and its first two derivatives continuous across 
the mixed layer. The exact form of the profile is unim- 
portant. For a plot of the profile m = 0.9, see Fig. [| 

The basic flow U(y) is obtained by requiring the ve- 
locity and all relevant derivatives to be continuous at the 
edges of the mixed layer: 

U{y) = 1 - Gy 2 /2 , for y<p, (6) 

U(y) = U(p)-G I dyy/fi, for p<y< P + q,(7) 
Jp 

U(y) = G (1 - y 2 ) /2m, for y>p + q. (8) 
Here G is the streamwise pressure gradient. 



E p (x,y,t) = - [u p (x,y,t) 2 +v p (x,y,t) 2 ] . (10) 



We can express the mean (over x) density of the kinetic 
energy as follows: 

E p {y,t) = ( K E p (x,y,t)^ = E p (y)exp(2 7p t), (11) 

£p(v) = l(K(y)\ 2 + \v P (y)\ 2 ) . 

The physics of our phenomenon will be discussed in 
terms of the balance equation for the averaged distur- 
bance kinetic energy. Starting from the linearized Navier- 
Stokes equations for u p and v pi dotting it with the dis- 
turbance velocity vector, averaging over one cycle in x 
and using Eqs. (p|)-(|Il|) leads to 



2 lp £ p {y) = V • J p (y) + W p+ (y) - W p ^(y) , (12) 

where the energy flux J p (y) in the y direction, rates of 
energy production (energy taken up by the primary insta- 
bility from the mean flow) W p + (y) and energy dissipation 
(by the viscosity) W p - (y) are given by 

\uJy)p*Jy) + c.c.l 1 
J p {y) ee 1 pKSnP ^J L + ±^ y) V£ p (y), (13) 

W P+ (y) = -\u'(y)[u p {y)v;{ y )+c.c] , (14) 

W p -(y) ee fM^2k%{y) + l[\u' p (y)\ 2 + \v' p (y)\ 2 ] 

The superscript * denotes complex conjugate. To plot 
these functions we need to solve Eq. (||) as an eigenvalue 
problem, to obtain c p , j p , and 4> p (y) at given Re and k p . 
The value of c p determines the position of the critical 
layer. It is convenient to compute and compare the space 
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averaged production and dissipation terms T p+ and T p _ 
defined by: 

T p± = J W p± (y)dy I J £ p (y)dy . (15) 

The local production of energy can be positive or neg- 
ative, indicative of energy transfer from the mean flow 
to the primary disturbance and vice-versa respectively. 
The production in one region (where W p+ (y) > 0) can 
be partly canceled out by a "counter-production" in other 
region (where W p+ {y) < 0). 

The use of these measures can be exemplified with the 
neat fluid (m = 1.0 here). The laminar flow displays its 
first linear instability at a threshold Reynolds number 
of Re t h = 5772, which means that the total production 
T p+ across the layer becomes equal to the total dissipa- 
tion T p _ at this value of Re. Examining Fig. |J we can 
see that the disturbance kinetic energy is produced pre- 
dominantly within the critical layer, where the basic flow 
velocity is close to the phase speed of the disturbance, 
while most of the dissipation is in the wall layer. The 
balance is not changed significantly when the viscosity 
ratio is changed to 0.9 so long as the mixed layer is not 
close to the critical layer. There is a small region of pro- 
duction and one of counter-production within the mixed 
layer, whose effects cancel out, leaving the system close 
to marginal stability. 

We now turn our attention to Fig. |^, in which our main 
point is demonstrated. The Reynolds number is the same 
as before, but the mixed layer has been moved close to the 
critical layer. It is immediately obvious that the earlier 
balance is destroyed. The counter-production peak in the 
mixed layer is much larger than before, making the flow 
more stable. The wavenumber used is that at which the 
flow is least stable for the given Reynolds number at this 
p. For m = 0.9, the threshold Reynolds number is 46400. 
Fig. ^ shows the energy balances at marginal stability - 
the picture is qualitatively the same here as at Re » 5772 
for the neat fluid. 



B. The mechanism of stabilization 

The main factor determining the instability is the 
energy intake from the mean flow, which is driven by 
the phase change caused by the viscosity stratification. 
The dissipation on the other hand depends only on 
Reynolds number and does not respond disproportion- 
ately to changes in viscosity. In neat fluids, the term 
containing U"{y) in (0) is always of higher order within 
the critical layer. However, with the introduction of a 
viscosity gradient within the critical layer, the gradients 
of the basic velocity profile will scale according to the 
mixed layer coordinate £. We show in the analysis that 
follows that for q < 0(Re -1 ^ 3 ), the term containing U" 
is now among the most dominant. Since most of the pro- 
duction of disturbance kinetic energy takes place within 
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FIG. 4: Energy balance: production W p +(y), solid line; dissi- 
pation W p -(y), dot-dashed line, Re = 5772. Top: m = 1, 
r p+ = T p - = 0.0148. Bottom: m = 0.9, p = 0.3, 
T p+ = 0.0158, T p _ = 0.0148. In this and the two subse- 
quent figures the solid vertical lines show the location y c of 
the critical lines, whereas the region between the dotted lines 
is the mixed layer. 
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FIG. 5: Energy balance: production W p +(y), solid line; dis- 
sipation W p -(y), dot-dashed line. Re = 5772, m — 0.9, 
p = 0.85, T p+ = -0.0114, r p _ = 0.0122. 



5 



o.io 



i 0.05 



0.00 



0.0 



0.2 



0.4 0.6 

v 



0.8 



FIG. 6: Energy balance: production W p +(y), solid line; dis- 
sipation W p -(y), dot-dashed line. Re = 46400, m = 0.9, 

p = 0.85, r p+ = rv = 0.0053. 



the critical layer, we return to equation (g) and isolate all 
lowest-order effects within the critical layer. The relevant 
normal coordinate in the critical layer is 



V 



y-y c 



(16) 



where y c is the critical point defined by U(y c ) = c, and e 
is the critical layer thickness, which is a small parameter 
at large Reynolds number. The basic channel flow veloc- 
ity may be expanded in the vicinity of the critical point 
as 



U(y) =c+(y- y c )U'(y c ) + {v ^- U"(y c ) 



(17) 



We use (|17|) , and redefine 4> p (y) = $(f?) and = u(£), 
to rewrite (|J) within the critical layer. We obtain 



e ~ Re~ 1/3 ee (fc p Re)~ 1/3 



(18) 



and the lowest order equation in the critical layer: 



IT) 



dy 



$"_-^ Xi /$ = v <S>^+2xv'<$>'"+x 2 v"<5>" , (19) 
v 



where \ = e /l is 0(\) for the mixed layer. In the absence 
of a viscosity gradient in the critical layer (i.e. v = 1), 
equation (fL9f) would reduce to 



IT] 



dU_ 
dy 



$" = $( 4 ) 



(20) 



which is the traditional lowest-order critical layer equa- 
tion for a parallel shear flow The mechanism for 
the stabilization now begins to be apparent: there are 
several new terms which can upset the traditional bal- 
ance between inertial and viscous forces. In order to 
narrow down the search further, we resort to numeri- 
cal experimentation, because although all terms in ( |T9| ) 



are estimated to be of O(l), their numerical contribu- 
tions are different. It transpires that the second term 
on the left hand side of ([l9]) is particularly responsible: 
it is straightforward to verify that it originates from the 
term containing U"(y) in the modified Orr-Sommerfeld 
equation. As testimony, note the dramatic effect on U" 
in Fig. 3. Any reasonable viscosity gradient of the right 
sign will pick up this term, leading to vastly enhanced 
stability. 

Indeed, in the light of this discussion we can expect 
that the large effect of retardation of the instability would 
even increase if we make the mixed layer thinner. This is 
indeed so. Nevertheless, one cannot conclude that insta- 
bility can be retarded at will, since other disturbances, 
differing from the primary mode, become unstable first, 
albeit at a much higher Reynolds number than the pri- 
mary mode; when we stabilize a given mode substantially, 
we should watch out for other pre-existing/newly desta- 
bilized modes which may now be the least stable. 

Finally, we connect our findings to the phenomenon of 
drag reduction in turbulent flows. Since the total dissi- 
pation can be computed just from the knowledge of the 
velocity profile at the walls, any amount of drag reduc- 
tion must be reflected by a corresponding reduction of 
the gradient at the walls. Concurrently, the energy in- 
take by the fluctuations from the mean flow should reduce 
as well. Indeed, the latter effect was measured in both 
experiments and simulations [[ll| The question 
is which is the chicken and which is the egg. In our cal- 
culation we identified that the reduction in production 
comes first. From Figs. 4 and 5 which are at the same 
value of Re we see that the dissipation does not change 
at all when the mixed layer moves, but the production is 
strongly affected. Of course, at steady state the velocity 
gradient at the wall must adjust as shown in Fig. 6. 



C. Continuous Viscosity Profile 

One could think that the strong stabilization discussed 
in the previous subsection is only due to the precise posi- 
tioning of the mixed layer at the critical layer. If so, the 
result would have very little generic consequence. In this 
subsection we show that any reasonable viscosity profile 
achieves the same effects. To this aim we consider the 
effective viscosity profile reported in Q (in their Fig. 5) 
which is obtained from simulations of a turbulent channel 
flow with polymer additive. It may be prescribed as 



Kv) = 1 + (m- 1) 



y-p 



(21) 
(22) 



with q ~ 0.4, and m ~ 0.7, as shown in Fig. [?]. The 
energy balance for the least stable primary mode at 
5772 for this case (Fig. H) shows a large counter- 



Re 



production of disturbance kinetic energy, which is in fact 
more pronounced than what we obtained with a mixed 
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FIG. 7: Prescribed continuous viscosity profile (in accordance 
with that obtained in direct numerical simulations of poly- 
meric flow). The corresponding laminar velocity profile U(y) 
and its second derivative are also shown. 



III. SECONDARY INSTABILITIES 

A laminar flow through a channel is linearly unstable at 
Re = 5772. In all except the cleanest experiments, how- 
ever, the flow becomes turbulent at much lower Reynolds 
numbers, as low as 1000 (IJ, This is because the 

linear stability analysis is carried out on a steady lam- 
inar velocity profile, whereas a real flow, except under 
carefully designed clean conditions, consists in addition 
of small but finite disturbances (most of whom will de- 
cay at long times). The stability behaviour of the real 
flow is quite different from that of the steady profile: 
the actual flow is unstable to new modes, often referred 
to as secondary modes. The secondary modes are often 
three dimensional, and their signature is prominent in 
fully-developed turbulence. As described below, the sec- 
ondary instabilities are studied by a Floquet analysis of 
the periodic primary flow we obtained earlier. 

As is usual in the analysis of secondary instabilities 
pHI |l7| ], we begin by splitting the flow into a periodic 
component (consisting of the mean laminar profile in ad- 
dition to the primary wave) and a secondary disturbance, 
e-g-, 
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FIG. 8: Energy balance: production W p +(y), solid line; dissi- 
pation W p -(y), dot-dashed line. Re = 5772, m = 0.7, p = 0.6, 
q = 0.4, T p+ = -0.0345, T p - = 0.0138. 



layer (Fig. ||). Thus the strong stabilization effect does 
not require careful placing of the viscosity variation at a 
particular layer. It is sufficient that there exist a viscos- 
ity variation in the region of the critical layer (indicated 
as the vertical line in Fig. ||) to achieve the stabilization. 

It comes as no surprise that this continuous viscosity 
profile behaves very similarly to the thin mixed-layer. If 
we return to equation we will see that all we have 
now done is to increase both v 1 (which is proportional to 
m — 1) and q threefold (the effective q here is closer to 0.3 
than 0.4, as we can see from Fig. |?]), so the ratio remains 
the same. 



where 



U(x,y,t) = U(y)x 

+A p (t) {[u p (y)x + v p (y)y] exp [ik p (x 



c p t)} 



(23) 



(24) 
c.c.} . 



Here x and y are units vectors in the x (steamwise) 
and y (wall normal) directions. The amplitude A p of 
the primary disturbance changes very slowly with time, 
and dAp/dt may be neglected during one time period. 
The spatial and temporal dependence of the secondary 
disturbance is written in the form 

u s (y, r±,t) = Ve^u s+ (y) exp [i (fe + • r± - uj+t)] 

+u s -(y)exp[i(k- -rj_ j, (25) 

where r± = xx + zz, and k± = k±x±k z z. We substitute 
the above ansatz into the Navier-Stokes and continuity 
equations, and retain linear terms in the secondary. On 
averaging over x, z and t, only the resonant modes sur- 
vive, which are related by 



fc 4 



therefore k± = ±q H — -x 



(26) 



for any vector q, and 



uj+ = uj + i7 s 



and 



w_ = (ujp + vy s . (27) 



Eliminating the disturbance pressure and streamwisc 
component of the velocity, we get the equations for the 
secondary disturbances v s and w s . Using the operator D 
for differentiation with respect to the normal coordinate 
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and the notation f± = —iw s ±/k z , the equations read 



0.02 



[i(w+ - fc+JJ) + i/(Z? 2 - fc+ - fc 2 ) + (£)!/)£)] 
x [(4 + fc 2 ) /+ - Dv+] - ik+U'v+ 

[ik + u p D + v p D 2 + ik_Dup\ v*_ 



(28) 




- k-k+) v p D + ifc+ (k 2 _ + k 2 z ) up] fl\=0, 



and 



[i(w+ - fc+17) + i/(Z? 2 



fc 2 - 



k\) + (D!/)D] (Df+ - v+) 



- ik+(DU) + {D 2 v)D + {Dv){D 2 - fc 2 - fc 2 ) 
Ap(kp + k—j 



iup (vt + Dfl) - -^Dv*_ 



Ap 
~2 



^ + D 2 ) i/, (IV,) 



S- =0, 



The boundary conditions are 







at 



y = ±1. 



(30) 



Equations (28) and (129), along with two corresponding 
equations in u* and fl, describe an eigenvalue prob- 
lem for the secondary instability. The four equations are 
solved by a Chebychev collocation spectral method, de- 
tails of the solution procedure are available in pl| . 

The most unstable secondary mode in our case is found 
to be the subharmonic, for which q = k z z. The produc- 
tion and dissipation are computed as before. 

We survey in turn the thin mixed-layer profile, and the 
continuous viscosity profile to see what viscosity variation 
does to the secondary instability. 



A. Mixed Layer 

The velocity and viscosity profiles here are as given in 
Fig. pi and the primary instability is that presented in 



Sect. 



Since the subharmonic (fc+ = fc- = k p /2) 



is the least stable mode, we present this case alone. In 
Fig. H a typical dependence of the growth rate of the sec- 
ondary mode on the spanwise wavenumber is shown. We 
can see that the viscosity variation damps the secondary 
mode significantly, but it is still unstable. However, there 
is a crucial difference in the primary instabilities of the 
two: the primary is unstable for a constant viscosity flow, 
but very stable in the mixed layer case. Therefore at long 
times, the secondary mode, which feeds on the primary 
for its existence, dies down in the latter case. To compute 
the time dependence of the amplitude of the secondary 
mode we computed the growth rate "f s by neglecting the 
time dependence of the amplitude of the primary words. 
As a result we obtain the growth rate 7 s L4 p (i)], in which 







/ \ 












/ \ 












/ / \ N 




/ / \ \ 




/ / \ \ 




. / \ \ 




'/ \ N - 




// \ N 




./ \ N - 




// \ \ 








// \. 











FIG. 9: Dependence of growth rate on spanwise wavenum- 
ber. Solid line: varying viscosity (p = 0.8, q = 0.1, m = 0.9); 
dashed line: constant viscosity (m = 1). k p — 1,A P — 
0.005, Re = 6000. 



(29) 



A p (t) can be an exponentially growing or a decaying func- 
tion of time. Having this growth rate we can present 
the time dependence of the amplitude of the secondary 
mode, see Fig. [lO]. Without the viscosity contrast, the 
amplitude of the secondary mode increases (essentially 
exponentially). With the viscosity contrast the ampli- 
tude decays in time. 

We now observe the balances of energy initially and at 
a later time in Figs, [ll] and |l2| respectively. The initial 
balance of energy is not very different from the constant 
viscosity case. At the later time, however, the production 
of secondary kinetic energy is significantly lower. The 
location y c of the critical point is seen from the figures 
to be close to the layer of stratified viscosity. If the two 
were well-separated, the stratification would do nothing 
to the secondary mode. 

A lowest-order analysis of the secondary stability equa- 
tions is not as straightforward as for the primary mode, 
since the secondary is highly dependent on the amplitude 
of the primary |lS| . We may however make the following 
observations from a critical layer analysis of equations 
( p8| ) and (^9|) and their counterparts. When A p 3> e, (cf. 
Eq.(|l8|)) only the nonlinear terms appear at the lowest or- 
der, and the secondary mode is completely driven by the 
primary. When A p ~ O(e), both the basic terms and the 
nonlinear terms contribute at the lowest order. It may 
be numerically determined, however, that the secondary 
is slaved to the primary here as well. When A p = o(e), 
the lowest-order theory for the secondary is (not surpris- 
ingly) exactly that given by ( |l^ ) for the primary. 

A direct estimate of the effect of the viscosity stratifica- 
tion on the secondary mode is obtained from the thresh- 
old amplitude A t h of the primary for the instability. At 
a Reynolds number of 6000 and primary wavenumber of 
k p = 1, for a neat fluid, all secondary modes are damped 
if A t h < 0.002, while for the continuous viscosity pro- 
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FIG. 10: Amplitude of the secondary mode in logarithmic 
scale as a function of time. Dashed line: constant viscosity, 
m — 1. Here 7 P = 0.0003, and the primary mode is unstable. 
Solid line: varying viscosity; here 7 P = —0.0206, the primary 
mode is stable. All conditions like in Fig. [| in particular 
A p (t = 0) = 0.005. 



FIG. 12: Production W a + and dissipation W s - of the kinetic 
energy of the secondary disturbance at time=40. Solid line: 
W a +,m = 0.9; dot-dashed line: W a -,rn = 0.9, A p = 0.00215; 
long dashes: W a +,m = 1; dotted line: W 3 -,m — 1,A P = 
0.00506. 
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FIG. 11: Production W a + and dissipation W s - of the kinetic 
energy of the secondary disturbance at time=0. Solid line: 
Ws+,m — 0.9; dot-dashed line: W a -,m = 0.9; long dashes: 
W a +,m — 1; dotted line: W a -,m — 1. The vertical lines 
show y c (the critical point location) for m = 0.9 (solid) and 
m — 1 (dotted). 



FIG. 13: Dependence of growth rate on spanwise wavenum- 
ber. Solid line: varying viscosity [according to equation (p^)]; 
dashed line: constant viscosity. Wavenumbers and Re as in 
Fig. | 



file, all secondary modes continue to be damped even for 
larger primary disturbances, up to A t h = 0.005. When 
the Reynolds number is reduced to 2000, the threshold 
amplitudes are 0.012 and 0.016 for the neat and viscosity- 
stratified fluids respectively. 



B. Continuous Viscosity Profile 

The velocity and viscosity profiles here are as given 
in Fig. [?], and the primary instability is that presented 



in Sect. II C. The counterparts for the continuous vis- 
cosity profile of Figs. [| to [l^ are presented in Figs. [l3| 
to |l6| respectively. It is clear that nothing has changed 
qualitatively. 

Fig. [jj] shows the dependence of the growth rate of the 
secondary mode on the amplitude of the primary distur- 
bance. It is clear that the instability is reduced by the 
stratification of viscosity, but there is no dramatic effect 
in the secondary alone. We may conclude that the large 
effect comes from the complete reliance of the secondary 
on the primary. 
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FIG. 16: Production and dissipation at time=40. Solid line: 
varying viscosity, A p — 0.00184; dashed line: constant viscos- 
ity, A p = 0.00506. 



FIG. 14: Variation ol the amplitude of the secondary in- 
stability mode with time. Solid line: varying viscosity, 
7 P = —0.0244; dashed line: constant viscosity, 7 P = 0.0003. 
Wavenumbers and Re as in Fig. Fol 
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FIG. 15: Production and dissipation at time=0. Solid line: 
varying viscosity; dashed line: constant viscosity, A v = 0.005 
for both. 



IV. CONCLUDING REMARKS 

We addressed the primary and secondary instability of 
simple channel flows, and examined the effects of small 
viscosity variations. We find dramatic effects of stabi- 
lization when the viscosity variations exist in the vicinity 
of the critical layers, in which the speed of propagation 
of the modes coincided with the mean velocity of the 
basic flow. With about 10% viscosity changes we can 



0.04 




FIG. 17: Dependence of the growth rate of the secondary 
mode on the amplitude of the primary disturbance. Solid 
line: stratified viscosity; dashed line: constant viscosity. k p — 
1, k+ = 0.5, k z = 1, Re = 6000. 



have very large increases in the threshold Reynolds num- 
bers for instability. In all cases we find that the main 
mechanism for the large effects is the reduction of the 
intake of energy from the mean flow to the putative un- 
stabale modes, which therefore become stable. For the 
same Reynolds numbers in Newtonian fluids there is no 
such mechanism for stabilization and these flows will be- 
come turbulent. We would like to propose that similar 
effects should be examined in the case of turbulent drag 
reduction by polymer additives. 

We recognize that in a turbulent flow there are many 
more modes that interact, but we propose that a simi- 
lar mechanism operates for each mode at its critical layer, 
where both elastic and viscous effects determine the mean 
flow. The advantage of the present calculation is that we 
can consider all the putative unstable modes, and con- 
clude that with a viscosity gradient similar to that seen in 
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polymeric turbulent flows the linear threshold Reth goes 
up five times (to 31000). We note in passing that this 
effect had not been put to an experimental test, and it 
would be exciting to have a confirmation of our predic- 
tions by future experiments. For actual turbulent flows 
we will need first to identify what are the main modes 
that interact between themselves and with the mean flow. 
A significant numerical effort is required, but appears 
worthwhile due to the importance of the phenomenon of 
drag reduction, and its relative lack of understanding. 

We have demonstrated that the exact form of the vis- 
cosity profile is immaterial; a continuous profile of vis- 
cosity in the critical region behaves exaclty like a thin 
mixed layer. We have shown that the secondary three 
dimensional modes of instability are "slaved" to the pri- 
mary linear mode of instability: the mechanism which 
stabilizes the primary mode indirectly ensures that the 
secondary is damped out quickly. 



Finally we note that a linear disturbance can rear its 
head either in the form of the fastest growing (or slowest 
decaying) mode as considered here; or in non-modal form 
with a transient growth followed by long-term decay jl9) . 
The former situation will correspond to relatively high 
Reynolds numbers, or cleaner set-ups. We expect similar 
conclusions in the latter situation as well. 
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